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Abstract 

We consider a two dimensional reflecting random walk on the nonnegative in- 
teger quadrant. This random walk is assumed to be skip free in the direction to 
the boundary of the quadrant, but may have unbounded jumps in the opposite 
direction, which are referred to as upward jumps. We are interested in the tail 
asymptotic behavior of its stationary distribution, provided it exists. Assuming the 
upward jump size distributions have light tails, we completely find the rough tail 
asymptotics of the marginal stationary distributions in all directions. This general- 



izes the corresponding results for the skip free reflecting random walk in 25]. We 



exemplify these results for a two node network with exogenous batch arrivals. 
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1 Introduction 



We consider a two dimensional reflecting random walk on the nonnegative integer quad- 
rant. This random walk is assumed to be skip free toward the boundary of the quadrant 
but may have unbounded jumps in the opposite direction, which we call upward jumps. 
Here, the boundary is composed of the origin and two half coordinate axes, which are 
called boundary faces. The transitions on each boundary face are assumed to be homo- 
geneous. This reflecting process is referred to as a double M/G/l-type process. This 
process naturally arises in queueing networks with two dimensional compound Poisson 
arrivals and exponentially distributed service times. Here, customers may simultaneously 
arrive in batch at different nodes. It also has own interest as a multidimensional reflecting 
random walk on the nonnegative integer quadrant. 
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A stationary distribution is one of the most important characteristics for this reflect- 
ing random walk in application. However, it is known very hard to analytically derive 
it except for some special cases. Thus, theoretical interest has been directed to its tail 
asymptotic behaviors. Borovkov and Mogul' skii [5] made great contributions to this prob- 
lem. They proposed the so called partially homogenous chain, which includes the present 
random walk as a special case, and studied the tail asymptotic behavior of its stationary 
distribution. Their results are very general, but have the following limitations. 



(a) The results are not very explicit. That is, it is hard to see how the modeling 
primitives, that is, the parameters which describe the model, influence the tail 
asymptotics. 

(b) The tail asymptotics are only obtained for small rectangles. No marginal distri- 
bution is considered. Furthermore, some extra technical conditions which seem to 
be unnecessary are assumed. 



For the skip free two-dimensional reflecting random walk, these issues have been ad- 



dressed by Miyazawa [25|, a nd its tail asymptotics have been recently answered for the 



marginal distributions in [20] , which considers the stationary equation using generating or 



moment generating functions and applies classical results of complex analysis (see, e.g. 



26| and references therein). This classical approach has been renewed in combination 



with other methods in recent years (e.g., see 19|, |24|). However, its application is limited 



to skip free processes because of technical reasons (see Remark 12 .31) . 

In this paper, our primary interest is to see how the tail asymptotics are changed when 
upward jumps may be unbounded. For this, we consider the asymptotics of the stationary 
probability of the tail set: 

{(it, i 2 ) e Z 2 + ; cth + c 2 i2 > x} 

as x goes to infinity for each C\,c 2 > 0, where Z + is the set of all nonnegative integers. In 
this paper, c = (ci,c 2 ) is called a directional vector if c\ + c\ = 1. We are interested in 
its decay rate, where a is said to be the decay rate of function p(x) if 

a — — lim — \ogp(x). 

n— >oo x 

Thus, we aim to derive the tail decay rates of the marginal distributions in all directions. 
These decay rates will be geometrically obtained from the curves which are produced by 
the modeling primitives, that is, one step transitions of the reflecting random walk. In 
this way, we answer our question, which simultaneously resolves the issues (a) and (b). 

Obviously, if the tail decay rates are positive, then the one step transitions of the 
random walk must have light tails, that is, they decay exponentially fast. Thus, we assume 
this light tail condition. Then, the mean drifts of the reflecting random walk in the interior 
and on the boundary faces are finite. Using these means, Fayolle, Iasnogorodski and 
Malyshev [12| characterize stability, that is, the existence of the stationary distribution. 
If the mean drifts vanish in all directions, then it is not hard to see that the stationary 
distribution does not have light tails in all directions. We also formally verify this fact. 
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Thus, we assume that not all the mean drifts vanish in addition to the light tail condition 
for the one step transitions. 

Under these conditions, we completely solve the decay rate problem for the marginal 
stationary distribution in each direction, provided the stability conditions hold. The decay 
rate may be considered as rough asymptotic, and we also study some finer asymptotics, 
called exact asymptotics. Here, the tail probability is said to have exact asymptotics for 
some function / if the ratio of the tail probability at level x to f(x) converges to a positive 
constant as x goes to infinity. In particular, if / is exponential (or geometric), it is said 
to be exactly exponential (or geometric). We derive some sufficient conditions for the tail 
asymptotic to be exactly exponential. 

The difficulty of the tail asymptotic problem mainly arises from reflections at the 
boundary of the quadrant. We have two major boundary faces corresponding to the coor- 
dinate axes. The reflections on these faces may or may not influence the tail asymptotics. 
Foley and McDonald 16| classified them as jitter, branch and cascade cases. We need to 



consider all of those influence to find the tail asymptotics. This problem becomes harder 
because of the unbounded jumps. 

To overcome these difficulties, we here take a new approach. Based on the stability 
conditions, we first derive the convergence domain of the moment generating function of 
the stationary distribution. For this, we use a stationary inequality, which was recently 
introduced by the second author 26] (see also @]), and a lower bound for the large 



deviations of the stationary distribution due to Borovkov and Mogul' skii [5j. Once the 
domain is obtained, it is expected that the decay rate would be obtained through the 
boundary of the domain. However, this is not immediate. We need sharper lower bounds 
for the large deviations in coordinate directions. For this, we use a method based on 
Markov additive processes (e.g., see (28|). 



We apply one of our main results, Theorem 13. 2\ to a batch arrival network with two 
nodes to see how the modeling primitives influence the tail asymptotics. We also show 
that the stochastic upper bound for the stationary distribution obtained by Miyazawa 



and Taylor [27| is not tight unless one of the nodes has no batch arrival. 



This paper is made up by six sections. In Section[2], we formally introduce the reflecting 
random walk, and discuss its basic properties including stability and stationary equations. 
In Section [3], main results on the domain and tail asymptotics, Theorems 13 . 1\ I3~2l and 13.31 
are presented. We also discuss about linear convex order to see the influence of the 
modeling primitives to the decay rate. Those theorems are proved in Section HI We apply 
them to a two node network with exogenous batch arrivals in Section [51 We give some 
remarks on extensions and further work in Section [HI 



2 Double M/G/l-type process 

Denote the state space by S = 7j 2 + . Recall that Z + is the set of all nonnegative integers. 
Similarly, Z denotes the set of all integers. Define the boundary faces of S as 

So = {(0,0)}, S 1 = {(i,0)eZ 2 + ;i>l}, S 2 = {(0,i)eZ 2 + ;i>l}. 
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Let dS = Uf =0 Si and S+ = S \ dS. We refer to dS and S+ as the boundary and interior 
of S, respectively. 

Let {Y(£); £ = 0, 1, . . .} be the random walk on Z 2 . That is, its one step increments 
Y{£ + 1) — Y(£) are independent and identically distributed. We denote a random vec- 
tor subject to this distribution by X^ + \ We generate a reflecting process {L(£)} = 
{(Li(^), L 2 (£))} from this random walk {Y(£)} in such a way that it is a discrete Markov 
chain with state space S and the transition probabilities p(i,j) given by 



P{L{£+l)=j\L{£) = i) 



P(X^=j-i), jeS,ieS + , 
P(XW=j-i), j E S,i E Sk,k = 0, 1, 2, 



where X^ is a random vector taking values in Z 2 . For this process to be non defective, 
it is assumed that X (+) > (-1,-1), > (-1,0), X (2) > (0,-1), and X (0) > (0,0) 
respectively. Here, inequalities of vectors are meant to those in component-wise. 

Thus, L(£) is reflected at the boundary dS and skip free in the direction to dS. 
In particular, their entries Li(£) and L 2 (£) have similar transitions to the queue length 
process of the M/G/l queue. So, we refer to {L(£)} as a double M/G/l-type process. 

Let R be the set of all real numbers. Similar to Z + , let M + be the all nonnegative real 
numbers. We denote the moment generating functions of X^ and X^ by 7 and 7^, 
that is, for 6 = (6*i, 62) G M 2 , 
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'0 ) = E {e^ x{+ \ lk {0) = E{e^ x(k) )), k = 0,1,2, 



where (a, b) is inner product of vectors a and 6. As usual, M 2 is considered to be a metric 
space with Euclidean norm ||a|| = ■sj (a, a). In this paper, we assume that 

(i) The random walk {^K(£)} is irreducible and aperiodic. 

(ii) The reflecting process {L(£)} is irreducible and aperiodic. 

(hi) For each e R 2 satisfying B\ > or 9 2 > 0, there exist t > and > such that 
1 < 7(£0) < 00 and 1 < 7&(^0) < 00 for = 0, 1, 2. 

(iv) Either £(xJ +) ) ^ or £(A^ +) ) ^ for = (x( +) , X { 2 +) ). 



The conditions (i) and (iii) are stronger than what are actually required, but we use them 
for simplifying arguments. For example, except for the exact asymptotics, (iii) can be 
weaken to the following condition (see Remark 12. 2p . 



(hi)' 7(0) and 7^(0) for k = 0, 1, 2 are finite for some 6 > 0. 



In the rest of this section, we discuss three basic topics. First, we consider necessary 
and sufficient conditions for stability of the reflecting random walk {L(£)}, that is, the 
existence of its stationary distribution, and explain why (iv) is assumed. We then formally 
define rough and exact asymptotics. 
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2.1 Stability condition and tail asymptotics 



It is claimed in the book of Fayolle, Malyshev and Menshikov 12] that necessary and 
sufficient conditions for stability are obtained (see Theorem 3.3.1 of the book). However, 
the proof of their Theorem 3.3.1 is incomplete because important steps are omitted. A 
complete proof can be found in 
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Furthermore, some exceptional cases are missing 
as we will see. Nevertheless, those necessary and sufficient conditions can be valid under 
minor amendment. 



In [12j, those necessary and sufficient conditions are separately considered according 



to whether all the mean drifts are null, that is, E(X\ 



E(X. 



or not. One may 



easily guess that the null drift implies that the stationary distribution has a heavy tail in 
all directions, that is, the tail decay rates in all directions vanish. We formally prove this 
fact in Remark 14.11 Thus, we can assume (iv) to study the light tail asymptotics. 



We now present the stability conditions of [12] under the assumption (iv). We will 
consider their geometric interpretations. For this, we introduce some notation. Let, for 
k 



1,2, 



m k = E(X 



k i 



m 



(i) 



E(X 



ft 



m 



(2) 



E(X, 



<2h 
k i 



Define vectors 



m 



(mi,m 2 ), 



m 



(i) 



(m 



(!) ™(1) 



1 5 



in, 



(i) 



i (i) w\ 



in, 



(2) 



— m 



(2) 
2 i 



m 

(2) 



(2) 



(2) „(2) 



1 J 



m\ ). 



Note that m ^ by condition (iv). Obviously, rn ® is orthogonal to for each 

k = 1,2. We present the stability conditions of 12| using these vectors below, in which 
we make some minor corrections for the missing cases. 



Lemma 2.1 (Corrected Theorem 3.3.1 of [12]) If m / 0, the reflecting random 
walk {Z(£)} has the stationary distribution, that is, it is stable, if and only if either 
one of the following three conditions hold. 



(I) mi < 0, wi2 < 0, (m, m-j 1 ^) < 0, and (m, m^'} < 0. 



,(2)\ 



(II) mi > 0, m 2 < 0, (m,m| L L ' ) ) < 0. In addition to these conditions, < is 

(2) 

required if m\ = 0. 



(2) 



in) nil < 0,m,2 > 0, (m,m±) < 0. In addition to these conditions, m^' < is 
required if m 2 = 0. 



(i) 



(2) (2) 

Remark 2.1 The additional conditions m\ < for m\ 



2 w „„ - in (II) and < for 

2 = in (El) are missing in Theorem 3.3.1 of 12]. To see their necessity, let us assume 

(2) 

m\ = in (II). This implies that the reflecting random walk can not get out from the 

(2) 

2nd coordinate axis except for the origin once it hits the axis. On the other hand, m\ 
may take any value because of no constrain on i n the first three conditions of (II). 



5 



Hence, 777.3 < is necessary for the stability. If = 0, then a similar situation occurs, 
but the first three conditions imply mf" < 0, and therefore no extra condition is needed 
for this case. By symmetry, the extra condition is required for (HI). It is not hard to see 



the sufficiency of these conditions if the proof in [12j is traced, so we omit its verification 



We will see that the stability conditions are closely related to the curves 7(0) = 1 and 
7fc(0) = 1. So, we introduce the following notation. 

r = {6 G R 2 ; 7(f) < 1}, dT = {6e R 2 ; 7 (0) = 1}, 

T k = {8e R 2 ; 7fc (0) < 1}, dT k = {0 E R 2 ; lk (6) = 1}, k = 1, 2. 



Remark 2.2 (a) V and F k are convex sets because 7 and 7^ are convex functions. Fur- 
thermore, T is a bounded set by the condition (i). (b) By the condition (iii), dT and 
dT k are closed curves. If we replace (iii) by the weaker assumption (iii)', then dT and 
dr k may be empty sets. In these cases, we redefine them as the boundaries of T and T^, 
respectively. This does not change our arguments as long as the solutions of 7(0) = 1 or 
7fc(0) = 1 ar e not analytically used. Thus, we will see that the rough asymptotics are still 
valid under (iii)', but this is not the case for the exact asymptotics. 



0-2 




Figure 1: Vectors for conditions (I) and (II) 

Note that m and are normal to the contact lines of dT and dT k at the origin 
and toward the outside of T and T k , respectively (see Figured]). From this observation, 
we can get the following lemma, which gives a geometric interpretation of the stability 
condition. We prove it in Appendix [A] 

Lemma 2.2 Under conditions (iii) and (iv), the reflecting random walk {L(£)} has the 
stationary distribution if and only if T PI T k contains a vector 6 such that 9 k > for each 
k — 1, 2. Furthermore, if this is the case, at least for either one of k — 1, 2, there exists a 
e r n T k such that 6 k > and 6 3 _ k < 0. 

Throughout the paper, we assume stability, that is, either one of the stability condi- 
tions (I), (II) and (HI), in addition to conditions (i)-(iv). We are now ready to formally 
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introduce tail asymptotics for the stationary distribution of the double M/G/l-type pro- 
cess. We denote this stationary distribution by v. Let L = L 2 ) be the random vector 
subject to the distribution u, that is, 

u(i) = P(L = i), ieS 

We are interested in the tail asymptotic behavior of this distribution. For this, we define 
the rough and exact asymptotics. We refer vector c G I 2 as a direction vector if ||c|| = 1. 
For an arbitrary direction vector c > 0, we define a c as 

a c = — lim — log P( (c, L) > x) , 

x^-oo x 

as long as it exists. This a c is referred to as a decay rate in the direction c. Thus, if 
the decay rate a c exists, P((c,L) > x) is asymptotically lower and upper bounded by 
e -(a c +e)x an( j e -(a c -e)x^ reS p ec tively, for any e > and sufficiently large x G K. If there 
exists a function / and a positive constant b such that P{{c, L) > x) ~ bf(x), that is, 

■im ^f.^U , (2.1) 

then P((c, L) > x) is said to have exact asymptotic f(x). In particular, if / is exponential, 
that is, f(x) = e~ ax for some a > 0, then it is called an exactly exponential asymptotics. 
It is notable that random variable (c, L) only takes a countable number of real values. 
Hence, we must be careful about their periodicity. 

Definition 2.1 A countable set A of real numbers is said to be 5-arithmetic at infinity 
for 5 > if 5 is the smallest number such that, for some Xq > 0, {x G A; x > x$} is 
a subset of {5n; n G Z+}. On the other hand, A is said to be asymptotically dense at 
infinity if there is a positive number a for each e > such that, for all x > a, \x — y\ < e 
for some y G A. 

The following fact is an easy consequence of Lemma 2 and Corollary in Section V.4a 



of 13|], but we prove it in Appendix [B] for completeness. 



Lemma 2.3 For a directional vector c > 0, let K c = {(c,n);n G Z^_}. Then, K c 
is asymptotically dense at infinity if and only if neither c\ nor c 2 vanishes and cyjc^ is 
irrational. Otherwise, K c is arithmetic at infinity. 

Because of this lemma, the x in (12 .ip runs over either arithmetic numbers 5n for some 
5 > or real numbers. Particularly, if K c is 1-arithmetic at infinity, then we replace x 
by n. For example, this n is used for the asymptotics: P{L^ = n,L 3 _k — i) ~ f(n,i) for 
each fixed % G Z+. 
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2.2 Moment generating function and stationary equation 



There are two typical ways to represent the stationary distribution of the double M/G/l- 
type process. One is a traditional expression using either generating or moment generating 
functions. Another is a matrix analytic expression viewing one of coordinates as a back- 
ground state. In this paper, we will use both of then because they have their own merits. 
We first consider the stationary equation using moment generating functions. Since the 
states are vectors of nonnegative integers, it may be questioned why generating functions 
are not used. This question will be answered at the end of this section. 

We denote the moment generating function of the stationary random vector L by 

ip(6)=E(e {e ' L) ), OeR 2 . 
We define a light tail for the stationary distribution v according to 26 . 



Definition 2.2 The v is said to have a light tail in all directions if there is a positive 
G R 2 such that f{0) < oo. Otherwise, it is said to have a heavy tail in some direction. 

Define the convergence domain T> of the moment generating function if as 

V = the interior of {0 G R; <p(6) < oo}. 

Then, we can expect that the tail asymptotic of the stationary distribution is obtained 
through the boundary of the domain T>. Obviously, D is a convex set because tp is a 
convex function on R 2 . Let us derive the stationary equation for this T>. Let 

<p+(0) = E(e< 9 > L h(L>0)) t 

<p k {6 k ) = E{e 6kLk l(L k > 1, L 3 _ fe = 0)), k = 1, 2, 

where l(-) is the indicator function. Then, 

<p{6) = <p + (0) + ipx{9 x ) + + <A)(0), 

where <£>q(0) = P(L = 0). From this relation and the stationary equation: 

L ~ L + X( + h(L e S+) + X^ k h(LeS k ), (2.2) 

fce{o,i,2} 

where ~ stands for the equality in distribution, and the random vectors in the right hand 
side are assumed to be independent, we have 

(1 - 7 (0)V+(0) = E - VMh) + (7o(0) - lVo(0), (2.3) 

fee{i,2} 

as long as <p(0) < 00. This equation holds at least for 6 < 0. 

Equation (12. 3p is equivalent to the stationary equation of the Markov chain {L(£)}, 
and therefore characterizes the stationary distribution. Thus, the stationary distribution 
can be obtained if we can solve (12.31) for unknown function <p, equivalently, <£>+, (fx, <p2 
and ipo. However, this is known as a notoriously hard problem. This is relatively relaxed 
when jumps are skip free (see, e.g., H, 2f|). This point is detailed below. 
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Remark 2.3 (Kernel method and generating function) To get useful information 
from (12. 3p . it is a key idea to consider it on the surface obtained from 1 — 7(0) = 0. 
This enables us to express <fi(9i) in terms of the other <£>-,( 6* j)'s under the constraint that 
1 — 7(0) = 0. Then, we may apply analytic extensions for (pi(9i) u sing complex variables 



for 0. This analytic approach is called a kernel method (see, e.g., [l9|, |24], |26| 



In the kernel method, generating function is more convenient than moment generating 
function. Let 7(21,22) = 7(l°g2i,log22), then 7(21,22) is the generating function cor- 
responding to 7(0). Note that 21227(21,22) is a polynomial of 21 and 2 2 . Particularly 
for the skip free two-dimensional reflecting random walk, 2i2 2 (l — 7(21,22)) = 0, which 
corresponds with 1 —7(0) = 0, is a quadratic equation of 2j for each fixed 2 3 _j. Hence, we 



can algebraically solve it, which is thoroughly studied in the book of Fayolle et al. [11 
However, the problem is getting hard if the random walk is not skip free. If the jumps are 
unbounded, the equation 1 — 7(21, 22) = has infinitely many solutions in the complex 
number field, and there may be no hope to solve the equation. Even if these roots are 
found, the analytic extension of fi(zi) would be harder. □ 



This remark suggests that the kernel method based on complex analysis is hard to use 
for the M/G/l-type process. We look at the problem in a different way, and consider the 
equation 1 — 7(0) = in the real number field. In this case, it has at most two solutions 
of 9i for each fixed 9 3 _i because 7(0) is a two variable convex function. However, we have 
a problem on the stationary equation (12.31) because we only know its validity for < 0. 
To overcome this difficulty, we will introduce a new tool, called stationary inequality, and 
work on T and Tk for k = 1,2. For this, moment generating functions are more convenient 
because they make T and to be convex. This convexity may not be true for generating 
functions because two variable generating functions may not be convex. 

Although we mainly use moment generating functions, we do not exclude to use gen- 
erating functions. In fact, they are convenient to consider the tail asymptotics in coordi- 
nate directions. For other directions, we again need moment generating function because 
C1L1 + C2L2 may not be periodic. Thus, we will use both of them. 



3 Convergence domain and main results 

The aim of this section is to present main results on the domain T> and the tail asymptotics. 
They will be proved in Section HI We first give a key tool for finding the domain V, which 
allows us to extend the valid region of (12. 3ft from {0 G M?;0 < 0}. 

Lemma 3.1 For 6 M 2 , (p(0) < 00 and (12. 3 p holds true if either one of the following 
conditions is satisfied. 

([2k.) e r and (p k (0) < 00 for k = 1, 2. 
([2b) e T n T 1 and (p 2 (9 2 ) < 00. 
fl2b) e r n r 2 and ifx(9x) < 00. 
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This lemma is a version of Lemma 6.1 of 26(, and proved in Appendix O The lemma 
suggests that r, T D Ti and T fl T 2 are important to find for <p(0) to be finite. They are 
not empty by Lemma [2T2l Obviously, these sets are convex sets. We will use the following 
extreme points of them. 

0(*,max) = arg max 7(^,02) = 1}, = arg min {9 k]1 (9 x ,9 2 ) = 1}, 

(01,02) (01,02) 

(fc ' c) = arg sup {9 k ; OeTf] T k } } {k ' e) = arg max {9 k ; 6 G dT D <9r fc }, 

(01,02) ( 6,1 ' 6, 2) 

where c and e in the superscripts stand for convergence parameter and edge, respectively 
(see Figure [2]). Their meanings will be clarified in the context of the Markov additive 
process {Z^\i)}. From these definitions, it is easy to see that, for k = 1, 2, 

{k ' e \ 7fc (0 (fc ' max) ) > 1, 



g(k,c) 



Q(k,ma,x) 



Ik 



<0(fc,max)\ < ^ 



Using these points, we classify their configurations into three categories as: 

(Dl) 0f' c) < dt' c) and 9 { 2 ' c) < d^ c \ (D2) (2 ' c) < (1 ' c) , (D3) (1 ' c) < (2 ' c) , 

where we exclude the case that 9^'^ > 9 ( [' c ^ and O^^ > 9^'^ because it is impossible 
(see Figures [2] and [3]) . These categories have been firstly introduced in [25| for the double 
QBD processes, and shown to be useful for the tail asymptotic problems. 
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Figure 2: Typical figures for (Dl) 
Using this classification, we define vector r = (ti,t 2 ) as 



{e[ hc \e^ c) ) if (di) holds, 

( iM'\0?' c) ) if (D2) holds, 
{0 { i' C \U8i ,c) )) if (D3) holds. 

where £ x and £ 2 are defined as 

l x {9 2 ) = max{0; j(9, 9 2 ) = 1}, £ 2 (0O = max{0; 7 (0 1? 0) = 1}. 

We will also use the following notation corresponding with them. 

^(9 2 ) = rmn{9 1 ; 1 (9 1 ,9 2 ) = 1}, ^(9,) = mm{9 2 ; 1 (9 1 ,9 2 ) = 1}. 



(3.1) 



Then, the domain T> is obtained as follows. 
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Figure 3: Typical figures for (D2) and (D3) 

Theorem 3.1 Under conditions (i)-(iv) and the stability condition, we have 

^ = {6>er max ;0<T}, (3.2) 
where r max = {0 e M 2 ; < 30' G T}. Furthermore, for k = 1, 2, 

r fc = sup{# fc > 0; ^fe(^fc) < 00}. (3.3) 



Remark 3.1 This result generalizes Theorem 3.1 of [25j in two respects. First, the 
domain is fully identified. Secondly, the skip free condition in 25] is relaxed. 



We next consider the decay rate of P((c, L) > x) as x goes to infinity. In some cases, 
we also derive its exact asymptotic. From the domain D obtained in Theorem 13 .![ we can 
expect that this decay rate is given by 



OL, 



sup{x >0;icG T>}. 



(3.4) 



However, it is not obvious to verify this equality. This is because the convergence param- 
eter of the moment generating function of a distribution on [0, 00) does not answer the 
tail decay rate of the distribution (see Appendix ID]) . 



We simplify notation a Bk as a k for ei = (1, 0) and e 2 = (0, 1). Note that 



a 1 



n, e 2 (-ri) > 0, 

Pi, Un) < 0. 



a 2 



72, §i(T 2 ) > 0, 

02, ei(-r 2 ) < 0. 



(3.5) 



where (3k is the positive solution x of 7(xefc) = 1 (see Figure 0] below). We now present 
asymptotic results in two theorems, which will be proved in the next section. 



Theorem 3.2 Under conditions (i)-(iv) and the stability condition, we have, 



lim - log P(L k > n, L 3 _ k = i) 

n— >oo fl 



z e z+, k = 1,2. 



(3.6) 



Furthermore, if r k 7^ 0<>> max ) anc j ^ Markov additive kernel {An ,n > —1} is 1- 
arithmetic, then we have the following asymptotics for some constant b k i > 0. 

lim inf e Tkn P(L k > n, L 3 _ fc = i) > b ki , ieZ+,k = l,2. (3.7) 
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Figure 4: Typical figures for «i = r x and «i = f3i (the colored areas are the domains) 

In particular, for (Dl) with k = 1,2, (D2) with k = 2 and (D3) with = 1, this is refined 

as 



lim e Tkn P(L k > n,L 3 _ k = i) = c k i 2 G Z + for some c k i > 0. 



(3.8) 



Remark 3.2 (a) In the case (Dl), if r k = ef'^ and 0( fc < max ) ^ (fc ' e) , then we can show 



lim e Tkn P(L k > n, L 3 _ k = i) = 0. 



(3.9) 



This is immediate from Remark 13.31 (b). We conjecture that (13. 9p is also true for r k = 

0(fc,m«) and ^(fc.max) = (fc,e)_ ^ when wg app j y ^ method to the gkip free cas6) 

^fc(z) with complex variable z (or the corresponding generating function) has a branch 
point at z = r k for r k = 0jp max ) under the case (Dl), which is a dominant singular point 
of (fk{z) (see 2(|). Thus, the rightmost point #( fc ' max ) corresponds with a branch point of 
the complex analysis. This is also the reason why (13. 9p holds because P(L k = n, = i) 
has the exact asymptotics of the form n c e~ Tkn with c : 



■|, — | as n — > oo. 



Theorem 3.3 Under the same conditions of Theorem 13.21 we have, for any directional 
vector c > 0, 



lim — log P((c,L) > x) — —a c 

i->oo x 



(3.10) 



where we recall that a c = sup{x > 0; xc G V}. Further assume that 7(a c c) = 1, 
7fc(a c c) 7^ 1 and a c Cfc ^ r k for fc = 1, 2. Then, we have, for some positive constant b c , 



lim e acX P((c,L) >x) = b c 



(3.11) 



where x runs over {5n;n G Z + } if (c, X^} is ^-arithmetic, while it runs over R + if 
is non-arithmetic, which is equivalent to C2 7^ and C1/C2 is not rational. 



Remark 3.3 (a) Since a k may be less than r k , the decay rate of P(L k > n) may be 
different from that of P(L k > n,L 3 _k = 0). (b) Similar to Remark 13.21 if a k = Q^ ,max > 
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and (fc ' max) ^ (fc ' e) , then we can show that <^(0 (fc ' max) ) < oo using the expression (jOg) 
in the proof of this theorem. This implies that 



lim e akn P{L k > n) 

n— >oo 



0. 



(3.12) 



This obviously implies f !3.9p . (c) If the jumps are skip free, then finer exact asymptotics are 
obtained for coordinate directions in 20] , which partially uses the kernel method discussed 
in Section 12.21 One may think to apply the same approach as in [20| . Unfortunately, this 
approach does not work because of the same reason discussed in Remark [ 



We now consider how the decay rates 77 and a c are influenced when the modeling 
primitives are changed. Obviously, if 7(0), 7i(#) and 72(0) are decreased by changing the 
modeling primitives, then the open sets T, Ti and T2 are diminished. Hence, we have 

Lemma 3.2 Under the assumptions of Theorem 13.21 if the distributions of X^ + \ X^ 
and X^ are changed to increase 7(0), 7i(0) and 72(0) for each fixed e R 2 , then the 
decay rates 77 and a c are decreased. 



Here, decreasing and increasing are used in the weaker sense. This convention will be 
used throughout the paper. To materialize the monotone property in Lemma 13.21 we use 
the following stochastic order for random vectors. 



Definition 3.1 (Miiller and Stoyan [30|) (a) For random variables X and Y, the dis- 
tribution of X is said to be less than the distribution of Y in convex order, which is denoted 
by X < cx Y, if, for any convex function / from R to R, 

E(f(X)) < E(f(Y)), 

as long as the expectations exist, (b) For two dimensional random vectors X and Y, if, 
(0, X) < cv (0, Y) for each e R 2 , then the distribution of X is said to be less than the 
distribution of Y in linear convex order, which is denoted by X <i cx Y. 



For linear convex order, Koshevoy and Mosler 22] gives several conditions to be equiv- 



alent (see also 34J|). Among them, the Strassen's characterization for convex order visu- 



alizes the variability of this order (see Theorem 2.6.6 of 30J). That is, X <i cx Y if and 



only if there is a random variable Ug for each e R 2 such that E(Ug\(0, X)) = and 

(0,Y) ~ (0,X) + U , (3.13) 

where ~ denotes the equality in distribution. 

Since e x is a convex function, the following fact is immediate from Theorems 13.21 and 
13.31 and Lemma [3.21 



Corollary 3.1 If the distributions of X w and X [2) are increased in linear convex 

order, then the decay rate and a c are decreased for k = 1, 2 and any direction c > 0, 
where the stability condition is unchanged due to (13.131) . 

This lemma shows how the decay rates are degraded by increasing the variability of 
the transition jumps X^ + \ X^ and X^K 
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4 Proofs of the theorems 



The aim of this section is to prove Theorems 13. 1[ 13.21 and 13.31 Before their proofs, we 
prepare two sets of auxiliary results. One is another representation of the stationary dis- 
tribution using a Markov additive process. Another is an iteration algorithm for deriving 
the convergence domain D. 



4.1 Occupation measure and Markov additive process 

We first represent the stationary distribution using a so called censoring. That is, the 
stationary probability of each state is computed by the expected number of visiting times 
to that state between the returning times to the boundary or its face. The set of these 
expected numbers is called an occupation measure. We formally introduce them. 

Let U be a subset of S, and let a u = mi{£ > 1;L(£) G U}. For this U, define the 
distribution g u of the first return state and the occupation measure as 

g u (i,j) = P (L{a u ) = j,a u < oo| L(0) = i) , i,j G U, 



a u -l 



h u (i,j) = E[ l(a u <oo) l^)=i) 



L{0)=i\, i,jeS\U. 



Let G u and H u be the matrices whose entries are g u (i,_i ) an d h u (i,j), respectively. 
H u may be also considered as a potential kernel (see, e.g., 32| for this kernel). Note 



that G u is stochastic since {L(£)} has the stationary distribution v. Let v v be the 
stationary distribution of G u , which is uniquely determined up to a constant multiplier 
by v v G v = u u . By E u u(a u ), we denote the expectation of a u with respect to u u . By the 
existence of the stationary distribution, E u u{a u ) is finite. Then, as discussed in Section 
2 of [29|], it follows from censoring on the set U that 

Kj) = ^Z7yE^« E P(i,k)h U (k,j), J£S\U } (4.1) 
vU ^ ' ieu k£S\u 

where p{i, k) = P (L(l) = k\ 2,(0) = i). 

We here need the distribution v v because U may not be singleton. The censored 
process is particularly useful when the occupation measure is simple. For example, if we 
choose dS for U, then h u (i, j) is obtained from the random walk {Y(£)}, which is simpler 
than {£(£)}. 

We next choose U — So U S2 for H u = {h u (i, j G S \ U}. This occupation 
measure has been used to find tail asymptotics of the marginal stationary distribution 



(see, e.g., |2l|, |25|, |28|). For each m, n > £ > 1, we denote the Z + x Z + matrix whose 
entry is h Ue ((m,i), (n,j)) by Hmn , where Ug = {0, 1, . . . , £ — 1} X Z+. Since this 
-^m(m+n) does not depend on m > £, we simply denote it by for n > 0, where Hq is 
the identity matrix. 
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We need further notation. For nonnegative integers m,n > 0, let 
r^((m, i), (m + n,j)) 

oo 

= Y / P(L(£) = (m + n,j), 
1=1 

m < Li(£) < min(L 1 (l),...,L 1 (£- 1))|X(0) = (m,i)). (4.2) 

For m > 1, rW((m, £), (m + n, j)) does not depend on m, so, for each n > 1, we denote 
the matrix whose (i, j) entry is r^((m, i), (m + n,j)) by i?„ for n > 1. On the other 
hand, for m = 0, we denote the corresponding matrix by for n > 1. 

For each n > 0, let fn be the vector whose i-ih entry is u(n,i). Similar to (14. ip . it 
follows from censoring with respect to U = U n that 

n-1 

Using the well known identity that 



n=0 \ n=l 

(14. 3 p can be written as 

n 

^ = ^Y.^^l t >u^ H m i} „> L (4 . 4) 



These formulas are also found in 28, 29]. Thus, Un can be also express by {R^} 
and {HrP}. These expressions will be also useful. 

The censoring process for U = So U S2 is a Markov additive process generated from the 
reflecting process {L(£)} by removing the boundary transitions on S U S , 2 . We will use 
this additive process in Section H~4"l So, we here introduce notations for it, and present a 
useful identity, called the Wiener-Hopf factorization. 

Denote this Markov additive process by {Z {1) (£)}. Specifically, let {A^;n > -1} be 
its Markov additive kernel, that is, entry of matrix An is defined as 

= P(Z«(£+ 1) = (m + n,j)\Z^\£) = (m,i)), n,m e G Z+, 

where the right-hand side is independent of m, then 



P(X (+) = (n,j -«)), i > 1, j > i — 1 
P(X« = (n,j-z)), z = 0,j>0 



1 

1 



for n > —1, where n and i, j are referred to as level and background states, respectively. 
Define matrix A* (i) for t > by 



n=-l 
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as long as they exist. This matrix is called a matrix generating function of the transition 
kernel Ail ■ Similarly, we denote the matrix generating function R£ by R^\t). Note 
that (gjZD is valid also for Z (1) (f) instead of L(£). Hence, R { *\t) is well defined for the 
Markov additive process {Z^\i)}. Then, we have the Wiener- Hopf factorization: 



I-A^(t)^(I-R^(t))(I-G^(t)), 



(4.5) 



where G*(t) is the Z x Z matrix whose entry is given by 

gS-\t)(ij) = e (t z ^°h(z?\a?) = j)\z£\o) = 

where cxf = inf{£ > 1; z[ (£) < 0}. The factorization (14.5j) goes back to [1], but t is 
limited to a complex number satisfying |t| < 1 for simplicity. The present version is valid 
as long as A^\t) exists. This fact is formally proved in 29], but has been often ignored 
(see, e.g., [28]) although it is crucial to consider the tail asymptotic 



in the literature 
problems. 



4.2 Iteration algorithm and bounds 

To prove Theorem 13.11 we prepare a series of lemmas. A main body of the proof will be 
given in the next section. We first construct a sequence in the closure of T>, denoted by 
V, and converges to r of (13. ip . For this, we rewrite (I2.3P as 

(l-7(%+(9) + (l-7k(%W 

= (73_ fe (0) - l)p 3 -k(03-fc) + (To(0) - %o(0). (4.6) 

An idea is to expand the region of <p(0) < oo using f l2.3l) and (14. 6ft with help of Lemma I3TT1 

Let rf ] = {0 G T k n T;6 3 - k < 0} for k = 1,2. Note that r^ 0) is not empty by 
Lemma \2. 21 Obviously, Lp^ k {9^ k ) is finite. Hence, by Lemma [3.11 <f+{0) and ip k (6 k ) are 
finite for G rf } . Thus, <p(0), </?i(0i) and (p 2 (6 2 ) are finite for G r[ 0) \JT ( 2 0) . We define 

0(A,O) _ ( ^(A,0)^(A,0) ) by 

e (A ' 0) = su V {9 k ;(9 l ,e 2 )eT^}, k = 1,2. 
By Lemma [2.21 at least one of 8[ A ' ^ and is positive under condition (iv). 

We next define {A ' n) = {d[ A ' n) , 6 (A ' n) ) for n > 1 by 

e[ A ' n) = su P {6 k ; eer k nr, e^ k < e^-\ 

Then, #( A ' ra ) is nondecreasing in n, and (A ' n) < max from our definition. Thus, the 
sequence #( A,ra ) converges to a finite positive vector because O^' 1 ' > 0. Denote this limit 
by 0( A '°°) = [o[ A '°°\ 6 , 2 A,0 °' ) ). Because T k (1 T is a bounded convex set, we can see that 

#( A <°°) = SU p{0 fc ; e g r fc n r, 3 - fe < 0£f >}> fc = 2 - ( 4 - 7 ) 

This can be considered as a fixed point equation, and we have the following solution, 
which is proved in Appendix [El 
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Lemma 4.1 O^'^ = r, and <p(0) is finite for all < (fc,oo) for k = 1,2, where 

0(i.°°) = (^ A ' oo) ,e 2 (^ A ' 00) )) ) (2,oo) = (^(^f' 00 ^,^* 00 ^- 

We need two more lemmas. We recall that c e I 2 is called a direction vector if 
||c|| = 1. Let 1 be the vector all of whose entries are unit. The dimension of 1 is either 2 
or oo, which can be distinguished in the context of its usage. 



Lemma 4.2 Let A(o) = {a + x G M 2 ; < x < 1} for a > 0. Then, under conditions 
(i), (ii) and (iii), we have, for any direction vector c > and any a > 0, 

liminf ~logP(L G xc+ A(o)) > - sup{(0, c); 7 (0) < 1}, (4.8) 

and therefore <£>(#) is infinite for G" r max , where r max is the closure of r max . 



Remark 4.1 This lemma is valid without condition (iv). Hence, it can be used for 
E(Xi) = E(X^) = 0. In this case, the right side of ( 14. 8 1) is zero, so the stationary 
distribution v cannot have a light tail. 

The first part of Lemma [4.21 is obtained in Theorem 3.1 of j^] (see also Theorem 1.6 
there). However, these theorems use Theorem 1.2 there, and its proof are largely omitted 
for the lower bound. Thus, the results are not well accessible, so we give its proof in 
Appendix [Fl We need one more lemma. 



Lemma 4.3 For each k = 1, 2, 



liminf — logP(Li > n, L% = i) > —6\ 

n— s-oo fl 



(l,c) 



i e Z.i 



(4.9) 



and therefore 9k > 0^'^ implies tpk(&k) = oo, and therefore y?(0) = oo. 

Proof. By symmetry, it is sufficient to prove the lemma for k = I. From (I4.4p . we have 



P(L 1 >n,L 2 = i ) = J2H ) ] l > 



(i) 



=71-1 



i e Z, 



Hence, (14. 9p follows from Theorem 4.1 of 21 



□ 



We are now ready to prove Theorem 13.11 



4.3 Proof of Theorem I3TT1 



We first prove that {0 G r max ; < r} C V. Since G r max implies the existence of 
0' G T> such that 0' > 0, this is sufficient to prove that 

V T = {0 G r ; < r} C V. 
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For cases (D2) and (D3), this V T is a subset of {6 G M 2 ; < O^'^} and {6 G M 2 ; 6 < 
O^ 1 ' 00 ^}, respectively, on which ip(0) is finite by Lemma 14.11 Hence, we only need to 
consider the case (Dl). In this case, tpk{9) < oo for 9 < 9^'°°^ = r k (k = 1,2), and 
therefore <p{0) < oo for 6 G V by Lemma [3.11 Thus, we have T> T C T>. Obviously, this 
also implies that 

r k = 9 { k k ' oo) < m V {9 k > 0; <p k (9 k ) < oo}. (4.10) 

We next prove that V C {8 G r max ; 6 < r}. Because of the symmetric roles of 9 1 and 
9 2 , it is sufficient to prove that either 9\ > T\ or 7(0) > 1 and g" r max implies <fi(9) = oo. 
The latter is immediate from Lemma 14.21 and therefore we only need to prove that 

Q\ > 7~i implies <pi(9i) = oo, (4-H) 

which together with (14. 101) also verify C I3 . 31) . We prove (14. lip for the cases (Dl), (D2) and 
(D3) separately. 

We first consider the case that (Dl) or (D3) holds. In this case, 

0(1,00) = e (l,c)^ 

therefore Lemma [4.31 verifies the claim. 

We next consider the case that (D2) holds. In this case, 9 ( {' ao ^ = £i(#2 2 ' c ^), and 9^^ = 
0(2,00)^ rp^g^ a pp}yi n g ^he same argument for 9i, we have <p{0) = oo and tp 2 {9 2 ) — 00 for 
9 2 > 9 { 2 ' oo) . If 9^' oo) = 9^' c \ then Lemma S3] again verifies the claim (14. lip . Hence, we 
assume that 9^'°°^ < 9± . 

In what follows, we consider the stationary equation ( 12.31) for 9\ G (0,9^'^). Let 
@2 = £ 2 (^i)- Since (f(0) < 00 for this 0, (12.31) is valid, and yields 

(i-7i(0i,£ 2 (0i))Vi(0i) 

= (72(^1, { 2 (^i)) - l)^(i 2 (9i)) + (70(01,^(01)) - l)Vo(O). (4.12) 

We increase 9i up to ^j 1 ' 00 ** in this equation. Note that 71 £ 2 (#i)) 7^ 1 and 72(^1, ^ 2 (^i)) 7^ 
1 for flf' min) < 9[ hc \ Let e = §(0?' c) - 0f ,oo) ). Since £ 2 "(#i 1,oo) ) = 0? ,oo) and~£ 2 (fl) is 
increasing for 9 G (6>[ 1 '°°' > , 6>| 1,c ^), we have, for any e G (0,e ), 

lim ^ 2 (£ 2 (#i)) = 00. 
This and (14.121) verify (14. lip . Hence, the proof of Theorem 13. II is completed. 

4.4 Proof of Theorem 13^1 

We prove Theorem 13.21 by two lemmas. They provide suitable bounds for the tail proba- 
bilities, and will immediately proves the theorem. 

Lemma 4.4 Under the assumptions of Theorem 13.11 

limsup — logP(Lfc > n, L 3 __fc = i) < — Tfc, i G Z+, k = 1, 2, (4.13) 
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and, for any directional vector c > 0, 

limsup -log P((c,L) >x)< - sup{u > 0; uc G V}. (4.14) 

Proof. (l4~T3j) is immediate from and (Q. To see f l4~T4j) . we use Markov 

inequality: 

e ua: P((c ) L) >x) < E(e {uc ' L} ), u > 0, n = 0, 1, . . . . 
Taking logarithm of both sides, dividing by x > and letting a; — > oo, we have (14.141) . □ 

Lemma 4.5 Assume the conditions of Theorem 13.21 Then, we have 

liminf - log P(Lx > n,L 2 = i) > -r ls i G Z+. (4.15) 

Further assume that the Markov additive kernel {A^; r7 > — 1} is 1-arithmetic concerning 
the additive component. If either (Dl) with t\ = 8^'^ < #( 1,max ) Q r (D3) holds, then, for 
some positive vector b, 

lim e Tin v^ ] = b, (4.16) 
while, if (D2) with t\ < O^'^ holds, then, for some constants b' > 0, 



liminf e rin i/ 1} > b'x, (4.17) 

n— >oo 

where the limit is taken component- wise and x is the left invariant vector of A* (ti). 



Proof. If n = 0^ 1,max ^ then we obviously have (14.151) by Lemma I4T21 Otherwise, (14.15[) 
follows from ( 14. 16ft and ( 14. 17ft or their 5-arithmetic versions for each positive integer 6. 
Thus, it remains to prove (I4.16P and (14.171) . 

We first consider the case where either (Dl) with 6^'^ < #( 1,max ) Q r (D3) holds. In this 
case, we have T\ = O^'^ = 6^' c ^ and therefore T\ < 6'j 1 ' max - ) by the definition ( 13. ip of T\. By 
Lemma 4.2 of 21], we already know that A^\e dl ) is positive recurrent. Furthermore, 
^2(62) < 00 for 6 2 < t 2 and 7i(0 (1 ' e ' ) ) = 1. Hence, using the notation Un of Section |2l we 
can verify all the conditions of Theorem 4.1 of [28] since is 1-arithmetic. Thus, we 

get Kim . 

We next consider the case where (D2) with n < 6?' c) holds. We use the same idea 
which is applied to the double QBD process, a special case of the double M/G/l-type 
process, in [2^ (see Proposition 3.1 there and Lemma 2.2.1 of j23[). However, the skip free 



condition for the reflecting process is crucial in the arguments of [23j, |25[ , and therefore 
we can not directly use it. Thus, we need some more ideas. 

Similar to the case (D3) for k = 1, (D2) implies that t 2 = and, for a positive 

(2) 

constant b\ for each I G Z + , 



lim e T2 V(£, i) = b 



(2) 
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Furthermore, let Xi be the z-th entry of the left invariant vector x of A* \e Tl ), where 



T2 = ^(ti), then it follows from Theorem A.l and (A. 11) of 2jJ that the complex variable 
generating function (z) = YltLo z ^ nas a simple pole at z = r 2 and no other pole on 
\z\ = t 2 because of the 1-arithmetic condition, and therefore, for some positive constant 

lim e T2i Xi = Cq. (4.18) 

i— >oo 

Hence, we have 

lim uPtyxr 1 = c lim e T2i v(£,i) = c bf\ £eZ + . (4.19) 

i— >oo i—toc 

Recall the Markov additive process {Z^\£)}, which is introduced in Section H~T1 We 
now change the measure of this process using T\ and x in such a way that the new kernel 
is defined as 

= A x \e^A^yA x , n>-l, (4.20) 

where is the diagonal matrix whose i-th entry is Xj. Obviously, An\i,j) is a non- 
defective Markov additive kernel. We denote the Markov additive process with this kernel 

by {Z (1 V)} = {{Z? (£),Z<P (£))}■ ^t 



i) 



n=-l 

which is the transition probability matrix of the background process {Z\ L> (n)}. By Lemma 



n=-l 

r(l)/ 



A. 2 of 21], A^ must be transient because T\ < 9^^ = 9i' e \ Under this change of 



measure, R { n ] and H { n ] are similarly changed to Rn and Hn\ Namely, 

^W = A-^e^W)^, n>l, ff« = K\^ n H^fA x n>0. 

It is notable that = Y^Li ^ s stochastic because x is also the left invariant 
vector of i?W(e Tl ) by the Wiener- Hopf factorization (14. 5p . For n = 1,2, . . ., let a^\o) = 
and let 

5 { l\n) = inf{£ > df\n - 1); z[ 1] (£) - z[ 1] (^\n - 1)) > 1}, 

then we can see from the version of (14.51) for the change of measure that Rn is the 
transition probability matrix at the first ascending ladder epoch, that is, its entry is 
given by 

rU(i,j) = P(zW(a! 1} (l)) - ^W(O) = n,zW(*P(l)) = j, \Z^(0) =i), n > 1, 

Because R^ is stochastic, this implies that Zi(a[ 1 \n)) drifts to oo as n — > oo with 
probability one. Let hn\i,j) be the entry of Hn \ Since 



hSKU) = l(n = 0)/ + ^^fW(z,A;)^W#, 
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we have, for n > 1, 

n 

= P(u? =1 {zf ViV)) - z{ 1} (0) = n} n {^ViV)) = j}I^ 1} (0) = 0, 

where the second equality is obtained because Z^\a^\tj) is increasing in t. Thus, 
hn is obtained as the solution of the Markov renewal equation, which is uniformly 

bounded by unit. However, we can not apply the standard Markov renewal theorem 
because its background kernel is transient. Nevertheless, we can show that, for some 
constant a > 0, 

lim Hjpl = -1. (4.21) 

n-^oo a 

Intuitively, this may be obvious because both entries of Z 1 (n) go to infinity as n — > oo 
and they behave like a random walk asymptotically when they get large. However, we 
need to prove (I4.2ip . Since this proof is quite technical, we defer it to Appendix iGl 

It follows from (I4.4p using the vector row y = {e~ T2i ; i > 0} that 

= [A x H^A x \u?yY . (4.22) 
Hence, applying the bounded convergence theorem and (I4.18p . we get f)4.17p . □ 



We are now ready to prove Theorems 13.21 and 13.31 

The proof of Theorem 13.21 

Because of symmetry, we only prove for k = 1. The rough asymptotic (13. 6p is imme- 
diate from Lemmas 14.41 and 14.51 The remaining part is also immediate from the second 
part of Lemma 14.51 □ 



4.5 Proof of Theorem 13^1 

We first consider where the ray xc with x > intersects T> in the (#i, #2)-pla n e- By 
Theorem I3.1[ there are three cases: 

(a) It intersects the vertical line Q\ — T\. 

(b) It intersects the horizontal line 82 = t~2- 

(c) It intersects dT (see Figure [5] below) . 
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(a) and (b) are symmetric, so we only need to consider cases (a) and (c). 
We first consider case (a). In this case, 



a, 



Cl 



-T\; 



because (c, c) = 1. On the other hand, it follows from ( 14. 15ft that 

lim inf — log P( (c, L) > x) > — lim inf — log P ( L± > — x, L<i = 

x->oo X C\ z-s-oo X V c l 

1 

> — n. 

Cl 



Hence, combining this with the upper bound ( I4.14p of Lemma [4.41 we have (I3.10p . 
We next consider case (c). In this case, (I3.10p is immediate from Lemmas 14.21 and 
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Figure 5: Typical figures for case (c) 



It remains to prove (13.111) . For this, we consider a singular point of the analytic 
function <f(zc) obtained from the moment generating function of (c, L). From (12.31) . 

2 

(1 - -f(zc))(p(zc) = 5^(7fc(zc) - -f(zc))ip k (c k z) + (7o(>c) - 7(zc))<p (0), (4.23) 

k=l 

for (3fiz)c G P. From the assumptions on a c , we first observe that the right hand side of 
(I4.23P is analytic for 3tz <= min(— Ti, ^2), which is greater than a c . 

We next consider the root of equation 1 — j(zc) = for $lz = a c . We claim that 
1 — j(zc) = has a single root at z = a c in the complex plane. From the above 
observation and the fact that j(zc) is an analytic function, the claim is proved. We 
consider the other root on = a c . Since j(0) is the moment generating function, we 
easily see that 

|7(*c)|< 7 ((9fk)c). 

Hence, if j(zc) = 1, then 7((9ffcz)c) > 1. The latter implies that ?R.z < or 3lz > a c since 
7(0) = 1 and j(xc) is a convex function for x G R. Thus, we must have dtz > a c for 
z E C satisfying 7(20) = 1. 
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Assume that (c, is (^-arithmetic for some 5 > 0. Since j(zc) is the moment gen- 

erating function of (c, X^), ^{zc) = 1 holds for Jfc > a c if and only if z = a c + 27ci/5, 
where i = y/—l. We change complex variable z to w = e 5z in y(^c), (fk(ckz), 'j(zc) and 
7fc(zc), and denote them, respectively, by f(w), fk(w), g(w) and gk(w), which are gener- 
ating functions. Since the analytic properties of the original functions are transformed to 
these generating functions, we can see that f(w) is analytic for \w\ < e Sac , and there is no 
singular point on the circle \w\ = e Sac except for w = e Sac . Since the circle in the complex 
plane is compact and g(w) = 1 has a simple root at w — e Sac , we can analytically expand 
f(w) given by 

f(w) = - l — — [ y2(g k (w) - g(w))f k (w) + (g (wc) - g(w))ip Q (0) ] 

to the region {w E C; \w\ < e s ° c + e} for some e > except for w = e Sac . Since h must 
be singular at w = e Soic and 1 — g(w) = has a single root there, f(w) has a simple pole 
at w = e Soic . Hence, we can apply the asymptotic inversion technique for a generating 
function of a 1-arithmetic distribution (e.g., see Theorem VI. 5 of [l4|). Thus, we have 

We finally assume that (c,X (+) ) is non-arithmetic. In this case, we back to moment 
generating function. Obviously, 3^2; = a c is impossible if z ^ a c because of the non- 
arithmetic condition. Furthermore, there is no sequence of complex numbers z n for n = 
1,2, .. . such that 7(z„c) = 1, lR,z n > a c and $lz n converges to a c as n — >• oo. This is 
because the existence of such a sequence implies that j(zc) = 1 for all zgC, which is a 
contradiction. This proves that 1 — 'y(zc) = has a single root at z = a c in the complex 
region {z G C; < < a c + e} for some e > 0. Hence, <f(zc) is analytically extendable 
for < ^tz < a c + e except a simple pole at z = a c . To this analytic function, we apply 
Lemma 3.2 of [3], which is an adaptation of the asymptotic inversion due to Doetsch [8|. 
Then, we can get (13.111) . 



5 Application to queueing networks 



In this section, we consider a two node Markovian network with batch arrivals. This 
network generalizes the Jackson network so that each node may have batch arrivals at 
once. By applying Theorem I3.2[ it is easy to compute the decay rates at least numerically 
not only for this modification but also for further modification such that the service rates 
are changed when either one of the nodes is empty because they can be formulated as the 
double M/G/l-type process. It may be notable that special cases of those models have 



been studied in the literature (e.g, see [l5|, |l7|, llSl, ll9| l2Q(). Exact asymptotics have been 



fully obtained there, while Theorem 13.21 only answers rough asymptotics. 

The aim of this section is twofold. First, we consider the influence of the variability of 
the batch size distributions to the decay rates. Secondly, we examine whether the upper 
bound of Miyazawa and Taylor [27| for the stationary distribution is tight concerning the 
decay rate, where it is assumed that there is no simultaneous arrival as in 27|. This 
tightness has been never considered. 
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5.1 Two node Markovian network with batch arrivals 



Consider a continuous time Markovian network with two nodes, numbered as 1 and 2. We 
assume that batches of customers arrive either at one of the two nodes or simultaneously 
at both nodes from outside according to the Poisson processes with rates A, and the 
service times at node k (k = 1,2) are independent and have the common exponential 
distribution with mean ^ , which are also independent of the batch arrivals. Because of 
this exponential assumption, service discipline is irrelevant as long as servers are busy. A 
customer who completes service at node 1 moves to node 2 with probability pi 2 , or leaves 
the network with 1 —p\ 2 . Similarly, a departing customer from node 2 goes to node 1 with 
probability p 2 i, or outside with probability 1 — p 2 \. To exclude trivial cases, we assume 



(Ek) < p 12 < 1 and < p 2 \ < 1. 



We assume without loss of generality that A + yUi + /i 2 = 1. 

By F, we denote the joint batch size distribution at each batch arrival instant. Let 
B = (B\, B 2 ) be a random vector subject to this distribution F, and denote its moment 
generating function by F. We assume 

([5b) For each non-zero 6 > 0, sup{t > 0; F{t6) < oo} = oo. 

Thus, F has light tails. The simplest model of this type is a two parallel queues with 
two simultaneous arrivals (no batch arrival at each node), whose exact tail asymptotics 



were studied by Flatto and Hahn |15|]. This result is recently generalized for a network 
but without batch arrival in [2(| . The present batch arrival network is more general than 
these model, but we can only derive the decay rates except for some cases. 

Let L t k be the queue length of node k at time t. Clearly, (L t i,L t2 ) is a continuous 
time Markov chain. We assume the intuitive stability condition: 

_ A(6i + 6 2 P2i) - _ X(b 2 + hp 12 ) 

(1 - Pl2P2ljA*l (1 - Pl2P2l)^2 

where bk = E(Bh). One can verify that this condition is identical with the stability 
condition given in Lemma 12.11 using fact that 

(m, mx ) = Ai2(A(6i + b 2 p 2i ) - - p 12 p 2 i)), 
(to, m { l ] ) = /Ui(A(& 2 + hp 12 ) - fi 2 (l - Pi 2 p 2 i)). 

Thus, (15.11) is indeed the stability condition for (L t i,L t2 ) to have the stationary distribu- 
tion, which is denoted by v. 

By the well known uniformization, we reformulate the continuous time Markov chain 
{{L t i, L t2 ); t > 0} as a discrete time one with the same stationary distribution v. This dis- 
crete time Markov chain is a double M/G/l-type process, so denoted by {(L^f), L 2 (£));£ = 
0,1,...}. 

This queueing network is more general than the model studied in [27j] in the sense 
that simultaneous arrivals at both nodes may occur. If there is no simultaneous arrival 
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at both nodes, then the model becomes a special case of the network of |27j because 
batch departures are not allowed. We will consider this special case in Section 15.31 for 



considering the quality of the upper bound of [27 



It is known that batch arrivals and/or simultaneous arrivals make it very hard to 
get the stationary distribution, while it is obvious if there is no such arrival. The latter 
network is the Jackson network, which has the stationary distribution of a product form 
as is well known. 



5.2 Influence of batch size distributions 

For the double M/G/l-type process for the batch arrival network, we compute 7(0) and 
7fe (0) for A; = 1,2. 

7 (0) = XF(9) + ^ ie - ei (i -p 12 + p 12 e 92 ) +/i 2 e- e2 (l - p 21 + p 21 e 9 ^), (5.2) 
7i(6>) = \F{0) + /i ie - ei (i - Pl2 + Pl2 e 92 ) + fx 2 , (5.3) 
l2 (0) = \F(0)+fi 1 +fi 2 e- 62 (l-p 21 +p 21 e e '). (5.4) 

Hence, if the distribution of B is increased in linear convex order, then 7(0), 71 (0) 
and 72(0) are increased for each fixed e R 2 as long as they exist. Hence, Lemma [3.21 
yields the following fact. 

Proposition 5.1 For the batch arrival Markovian network satisfying the conditions (Eh), 
([5b) and (15.11) if the distribution of B is increased in linear convex order, then the decay 
rates T\, t 2 in Theorem 13.21 and a c in Theorem 13.31 are decreased. 

To get the decay rates, we need to find (fc,e) and #( fc > max ) for k = 1,2, which are the 
roots of the equations 7(0) = 7fc(0) = 1 and 7(0) = 1 satisfying d g 9k - = 0, respectively. 
We here note that 7(0) = 71 (0) = 1 is equivalent to 71 (0) = 1 and 

e 92 = 1 -P2i+P2ie 9 \ 

which follows from 7(0) = 71 (0). Thus, the numerical values of the decay rates are easily 
computed using a software such as Mathematica, but their analytical expressions are very 
hard to get except for the skip free case. The model studied by Flatto and Hahn 15] is 
the simplest one among the latter case. Theorems 13.21 and 13.31 are fully compatible with 
their asymptotic results. 



5.3 Stochastic upper bound of Miyazawa and Taylor 

We next consider the stochastic upper bound for the stationary distribution u, obtained 



by Miyazawa and Taylor [27J. Since their model does not allow simultaneous arrival, we 
assume that either one of B\ and B 2 is zero. Thus, the joint batch size distribution F 
can be written as 

F(x u x 2 ) = F( Xl , 0) + F(0, x 2 ), x x ,x 2 > 0. 
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Let Fx(x) = F(x,0)/F(oo,0) and F 2 (x) = F(0, x)/F(0, oo), and let Ai = AF(oo,0) and 
A2 = AF(0,oo). For computational convenience, we switch to generating functions from 
moment generating functions. Let F k be the generating functions of Ff~. We present the 
upper bound of [27j using our notation. 

Proposition 5.2 (Corollary 3.2 and Theorem 4.1 of (27|) If flSa), fl5b) and the sta- 
bility condition (15 .ip hold, then the equations: 

Ai (Fi{si) - lj + /iis^^l - si) = H2P21S2 1 {1 - si), (5.5) 
A 2 (f 2 (s 2 ) - lj + H2S2 X (l - s 2 ) = /iiPi 2 s^ 1 (l - s 2 ), (5.6) 

have solutions (si, s 2 ) > 0. Let (/ii, h 2 ) be the maximal solution among them, then, 
P(L >n)< h^ ni h^ n \ n = (m, n 2 ) > 0, 

where L is a random vector subject to the stationary distribution v. 

To compare h k with the decay rate, we let 

r] k = \ogh k , fc = l,2. 

By Theorem 13.31 and Proposition 15.21 we have that % < a k . 

Letting s k = e dk in (15.21) . (15.31) and (15. 4p . 7(0), 7i(#) and 7 2 (0) can be written as 

AiFi(sx) +/iis^ 1 (l -p 12 +P12S2) + A 2 F 2 (s 2 ) + / u 2 s 2 1 (l -P21 +P21S1) = 1, (5.7) 
AiA(si) + /xisr^l - P12 + Pi 2 s 2 ) + A 2 F 2 (s 2 ) + M2 = 1, (5.8) 
AiFxCsx) +m + A 2 F 2 (s 2 ) +/i 2 s 2 - 1 (l -p 21 +P21S1) = 1- (5.9) 

Note that (15.51) and (15.61) imply (15. 7p . That is, (hi,h 2 ) satisfies the equation (15. 7p for 
variable (si, s 2 ). In other words, the point (771,77^ is on the curve <9r. 

Our question is when (r/i, r/ 2 ) is identical with («i, a 2 ), that is, when the upper bounds 
agree with the decay rates of the marginal distributions in coordinate directions. We 
answer it by the following theorem. 

Theorem 5.1 Under ([5b) and the stability condition (15.11) . the decay rate r\ k = 

logifc of the stochastic upper bound of [27] is identical with the decay rate a k for k = 1, 2 
if and only if both nodes have no batch arrival. 

Proof. The sufficiency of the single arrivals is immediate from the well known product 
form solution for the Jackson network. Thus, we only need to prove the necessity. Assume 
that (771 , 772) = (ocx,a 2 )- As we already observed, this implies that («i,a 2 ) G dT. Hence, 
from Theorem 13.31 and (13.51) . we can see that neither oti = 0i nor a 2 = /3 2 is possible 
because a k = (3 k implies that a 2 _fc = 0, where (3 k is defined after (13.51) . Furthermore, 
we cannot simultaneously have ot\ = #[ 1,max ) an d a 2 = 6^ 2 ' max ) because of («i,a 2 ) G dT. 
Hence, we must have either a\ = Q^^ or a 2 = 6* 2 2 ' e \ 
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Suppose that oci = 6^' e \ which implies that rji = O^'^ by our assumption. For 
convenience, we introduce notations: 



,(fc,e) e (k ' e) ■ , 10 

t\' J =e°' , i,k = 1,2. 



10ns 



Then, 771 = O^'^ is equivalent to hi = ti , and (t^ \ t^'^) is the solution of the equat 
(15. 7p and (I5.8P for variable (si, s 2 ). Thus, it follows from (15.71) and (15. 8p that 

^•^l-^l+^ltx 1 ^ (5.10) 

Substituting flBTTUD into (J521 with (s 1 ,s 2 ) = ($' e) , #' e) ), we have 

Ai (f^) - 1) + a 2 (f 2 (i -p 21 +p 21 4 1 ' c) ) - l) + /n(l -Pi^Kfif't 1 - 1) = 0. 

Thus, tj^ is obtained as the unique positive solution of this equation. Since hi = tf' c \ 
this implies that 

Ax (P^hx) - l) + A 2 (P 2 (l-p 21 +p 21 h 1 ) - l) +fi 1 (l-puP2i)(K 1 " 1) = 0. (5.11) 

From (15. 5p with (s x , s 2 ) = (/ii, /i 2 ) and (15. lip , we have 

A 2 F 2 (1 -p 2 i +P2ih) = nmmih^ 1 (1 - h{) - fx 2 p 21 h 2 \l - hi). 
This yields 

F 2 (l -P21 +P2ih) - 1 j j 

(1 - pal +P3lftl) - 1 

On the other hand, it follows from (15. 6 j) with (si, s 2 ) = (/ii, /i 2 ) that 

A 2 — — = A* 2 ft 2 - VxPvihi ■ 

h 2 - 1 

Hence, we must have 



F 2 (l - p 2 i + p 2 ih) - 1 = F 2 (h 2 ) - 1 

(1 - pii + i?2l/ll) - 1 

Since F 2 (s) is a strictly increasing convex function, this equation is true only when either 
F 2 (s) = s, which is equivalent to no batch arrivals at node 2, or 

h 2 = 1 -pn+pnhx. (5.12) 

Suppose that node 2 has batch arrivals. Then, (I5.12p holds, and therefore (I5.10p 
implies h 2 = 4 ■ This is equivalent to that (^1,^2) = d^- 1 '^. Hence, the assumption that 
(V11V2) = ( a i, a 2) implies that we must have the case (D2), which in turn implies that 
^ = 4^=4^). Hence, by the same arguments, we have either Fi(s) = s or 

hi = 1 -pi 2 + pi 2 h 2 . (5.13) 
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However, if both of ( I5.12p and fl 5 . 1 3 [) hold, then hi = h 2 = 1, which is impossible. Hence, 
we must have that Fi(s) = s, that is, node 1 has no batch arrivals. 

Applying Fi(hi) = hi to ( 15. 8ft with (si,s 2 ) = (hi,h 2 ) and using the fact that hi > 1, 
we have 

/iihi 1 = Ai + p 2 p 2 ih 2 x . (5-14) 
On the other hand, from ( 15. 6ft with (si, S2) = (/ii, ^2) and ( 15.121) . we have 

A 2 (F 2 (h 2 ) - l) = /iipia/ir^l " h 2 ) - ^h 2 \l - h 2 ) 
= (1 - /i 2 ) (nih^pu - ^hz 1 ) 
= (1 - h 2 ) ((Ai + li 2 Vi\h 2 v )vvi ~ ^h^) , 

where the last equality is obtained by substituting ( I5.14p . Rearranging terms in this 
equation and recalling that p 2 = (Ai + A 2 pi 2 )/((1 — Pi2P2i)A*i), we have 

A 2 {F 2 {h 2 ) - h 2 } = (1 - h 2 ) ((A 2 + A1P12) - (1 - ^21^12)^2^2 X ) 

= (1 - h 2 )(l - p 2 ip 12 )p, 2 (p 2 - h 2 l ) . (5.15) 

Since this p 2 is identical with the geometric decay rate of the Jackson network with single 
arrivals at both nodes, we obviously have p 2 < h^ 1 . Thus, the right-hand side of ( 15.151) 
is not positive, but its left-hand side must be positive because node 2 has batch arrivals, 
which implies F 2 (h 2 ) — h 2 > 0. This is a contradiction, and node 2 can not have batch 
arrivals. This concludes that both nodes do not have batch arrivals. By symmetry, we 
have the same conclusion for a 2 = 6^' e \ Hence, we have completed the proof. □ 



This theorem shows that the stochastic bound of 27j cannot be tight even for the 
decay rates. However, this may not exclude the case where one of the upper bounds 
is tight. In our proof, the tightness leads to a contradiction if either node 2 has batch 
arrivals or (D2) holds. This suggests that, if node 2 has no batch arrivals and if (D2) 
does not hold, then node 1 with rji = ai may have batch arrivals. The following corollary 
affirmatively answers it. 



Corollary 5.1 Under the same assumptions of Theorem 15.11 rji = ai holds if there is no 
batch arrival at node 2 and either (Dl) with ^( 1 ' max ) ^ f 1 or (D3) holds. These conditions 
are also necessary for rji = ol\ if node 1 has batch arrivals. Similarly, rj 2 = a 2 holds if 
and only if there is no batch arrival at node 1 and either (Dl) with #( 2 ' max ) g r 2 or (D2) 
holds. These conditions are also necessary if node 2 has batch arrivals. 



PROOF. By symmetry, we only need to prove the first two claims. If there is no batch 
arrival at node 2, then (hi, h 2 ) is obtained as the solution of the equations ( 15. 5ft and ( 15.61) 
with F 2 (s 2 ) = s 2 , that is, 

Ax (Fi{ Si ) - l) + fns^(l - hi) = p 2 p 2 is 2 \l - si), (5.16) 
P-2S2 1 = A 2 + P1P12S1 1 . (5.17) 
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Substituting /i 2 s 2 1 of (I5.17P into ( 15 . 1 6[) implies 

Ai (li(si) - l) - \2P21 (1 - s 1 ) + /n(l - pi 2 p 2 i)(l - Sl )sr a = 0. (5.18) 

On the other hand, this equation also follows from (15.71) . (I5.8P and the single arrivals at 
node 2, and therefore we have hi = ti\ equivalently, rji = 8^'^. Hence, rji = a\ holds if 
either (Dl) with #( 1 > max ) ^ p 1 or (D3) holds, because these conditions imply a.\ = 9^' e \ 
This proves the first claim. To prove the necessity, we note that T]i = a,\ implies a,\ = 9^' e \ 
Assume that node 1 has batch arrivals, then neither the batch arrivals at node 2 nor (D2) 
is possible as shown in the proof of Theorem 15. II Hence, it is required that node 2 has no 
batch arrivals and (D2) does not holds. The latter together with a% = 9^' e ' implies that 
either (Dl) with 0( 1 > max ) ^ f i or (D3) holds. Thus, the second claim is proved. □ 



6 Concluding remarks 

In this paper, we have studied the tail decay asymptotics of the marginal stationary dis- 
tributions for an arbitrary direction under conditions (i)-(iv) and the stability condition. 
Among these conditions, (i) is most restrictive for applications. For example, it excludes 
a priority queue with two classes of customers. However, it can be relaxed as remarked 
in Section 7 of j25|. Hence, (i) is not a crucial restriction. As we already noted, condi- 
tion (iii) can also be relaxed to (iii)' for obtaining the decay rate. Thus, the rough tail 
asymptotics can be obtained under the minimum requisites. 

What we have not studied in this paper is of other types of asymptotic behaviors of 
the stationary distribution. In particular, we have not fully studied exact asymptotics. 



We have recently studied this problem for the skip free reflecting random walk in [20 
For the unbounded jump case, this is a challenging problem. On the other hand, finer 
asymptotic results are available for some specific models in the literature (e.g., see Q), 
but require complex analysis. We hope the convergence domain which we obtained in this 
paper may be helpful for the kernel method as demonstrated in 2^ . 
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A Proof of Lemma 



2.2 



We first claim that, for each k = 1,2, there is a G T R such that > if and only if 

(0,m)<O, (6,m {k) )<0 and 4 > 0. (A.l) 

Because of symmetry, we only prove this for k — 1. To prove this, we define functions / 
and fi as 

/(«) = £(e"<^ (+) >), = E(e"<^ (1) >), u G K, 

as long as /(it) is finite for each fixed Obviously, by condition (iii), 

E{(6,XM)) = e 1 E(X { 1 +) ) + 9 2 E(X { 2 +) ) 

exists and finite for any G R. Choose a G R such that #i > and /(l),/i(l) < oo. 
This exists by (iii). Since /(0) = 1 and f(u) is convex in u, 

(e,m) = E((0,xM)) = f (0)<0 (A.2) 

is necessary and sufficient to have that /(ito) < 1 for some Uq > 0, which is equivalent to 
that uqO G T and u 9x > 0. Let 6' = u o 0, then (I A. 21) is equivalent to that (0 1 , m) < 0. 
Using this 6' , we apply the same arguments to function fx and can see that (6' , m^'j < 
holds if and only if there is a ux > such that U\& G r x and > 0. This implies that 
mm{\, U\)Q' G rnTi since T and Tx are convex sets. Thus, the claim is proved. 

We next show that (1A.1I) for k — 1 follows from either one of the stability conditions 
(I), (II) and (in). We first assume (I). Since > 0, we consider the possibility that 
= 0. In this case, m± < by the first inequality in (I) and m 2 < 0. Hence, 
Tx = {0 G R 2 ; <pi{0) < 1} is the region between two straight lines 9x = and 8x = a for 
some a > 0. Since T is not empty by (iii), we must have (IA.10 for some such that 9x > 
and 02 < 0. Thus, we can assume that rwp > 0. We put — {6x, 9 2) as 

{(!) (!) \ n 

-mi - e, > 0, 

— e, mj/ < 0, 

where e > is chosen so that (rn, mj 1 ) — em 2 < and m^mx — em 2 < 0, which is 
possible by (I) and mi < 0. Note that 9x > and #2 < in this definition. Then, we have 

(0, m) = m^mx - ((m^ + e)l(mS 1} > 0) + el(mf ) < 0))m 2 

= m.j^) — em 2 j l(m[ 1 ' ) > 0) + (m^Wi — em 2 )l(m^ < 0) < 0, 

(0, m (1) ) = m^m^ - ((m^ + e) 1(771^ > 0) + el(mi 1} < 0))m^ 1} 

= m^m^^m^ < 0) -emS 1 ' < 0. 

Thus, we have I 1A.1I) for = 1. 

We next assume (II). We consider the possibility that mi = 0. In this case, we put 
0=(9x,9 2 ) as 

X = -m 2 - e, # 2 = e, 
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where e > is chosen so that — m 2 — e > and (m, rrvp) — e(rrip — rnp) < 0, which is 
possible by (II). In this case, 6i,9 2 > 0. Then, we have 

(6, m) = — (m 2 + e)m\ + em 2 = em 2 < 0, 

(0, m^) = — (m 2 + e)m^ + em 2 = ( m i m ± ) — e{rrip — mp 1 ) < 0. 

Thus, we can assume that vii\ > 0. We choose e > such that (m,m ( f > ) + em\ < 0, 
which is possible by (II), and put 

9 1 = m£ ) + e, 6 2 = -m [ p. 

Then, 9\ > and 9 2 > since rrip > and mf < by (II). Hence, we have 

(6, m) = (m^ + e)mi — m^ m 2 = (m, Tnp) + e mi < 0, 
(0,m^) = {m 2 + e)mP — m^m^ = emf < 0. 

We finally assume (El). In this case, if m 2 = 0, then rap < 0, and therefore Ti is the 
region between #1 = and 9\ = b for some b > 0. Since mi < and m 2 > 0, we can easily 
see that (1A.1I) holds true for 6>i > and 9 2 < 0. Thus, we can assume that mP > 0, and 
therefore we can choose e > such that 

ernfp + Wbpm\ < 0, 

and put 6 = (e,mi). Then, 9\ > 0, 9 2 < 0, and 

(0, = e mi + m\m 2 < 0, 
(0, ) = e trip + mimj 1 ' < 0. 

Thus, we have shown (lA.ip for k = 1. Furthermore, 9 2 < for (I) and (HI). By symmetric 
arguments, (lA.ip for A; = 2 is obtained, and #i < for (I) and (II). Thus, either one of the 
stability conditions of Lemma 12.11 implies flA.ll) . The converse is immediate from flA.ll) 
and the observation that is presented just before this lemma. 



B Proof of Lemma 2.8 



Obviously, if either c\ or c 2 vanishes, then K c is arithmetic. Hence, we assume that c\ ^ 
and c 2 7^ 0. If ci/c 2 is rational, we obviously see that K c is arithmetic. Thus, we only 
need to prove that K c is asymptotically dense at infinity if ci/c 2 is irrational. For this, 
we combine the ideas which are used for proving Lemma 2 and Corollary in Section V.4a 



of 13 



Assume that Ci/c 2 is irrational and c\ < c 2 . The latter can be assumed without loss 
of generality because the roles of C\ and c 2 are symmetric. For each positive integer n, let 

A(n) = {cimi — c 2 m 2 ; < Cimi — c 2 m 2 < c 2 , m 2 < n, m 1 ,m 2 G Z + }. 
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then the number of elements of A(n) is strictly increased as n is increased because of the 
irrationality. Hence, for each e > 0, we can find positive integer n such that there are 
it, u' G A{n) such that \u — u'\ < e because A(n) is a subset of the interval [0,02]. Since 
we can find mi, rri2, m! 1 , m' 2 such that m\ > m^, u = c\m\ — C2rri2 and u' = C\vn! x — C2m' 2 , 
we have 

\c\{mi — m^) — c 2( m 2 — m 2) \ — \u — u'\ < e. 

Since nil > m^, we obviously require that m 2 > m 2 . Hence, we put a = 02(7772 — m' 2 ), 
then we have |x — y\ < e for some x, y G iY c and x,y > a. 



C The proof of Lemma 13.1 



We only prove this lemma when ([2^l) is satisfied because the other cases are similarly 
proved. We immediately have ( 12.31) if <p+(0) < oo. For proving this finiteness, we apply 
truncation arguments for (12. 2ft . For each 77 = 1, 2, . . ., let 

f n (x) = min(x,n), x G R. 

If x < 77, then / n (x + y) < x + 7/ = /„(x) + t/. Otherwise, if x > n, then f n (x + y) < n = 
f n (x). Hence, for any x > and y G R, 

/n(x + y)</„(x) + | £ (C.l) 

From (12. 2p . we have 

(0, j L)^(0, j L) + (0,X«>1( j LgS + ) + ]T (0,X«)1(XG^). (C.2) 

fce{o,i,2} 

Hence, we have, using the independence of L, X^°\ X^ and X^, 

E{e^°> L ») < E(e f "^h(L G 3+, (0, L) < n))E{e^ x(+) ») 
+E(e M{e ' L)) l(L G S+, (0, L) > 77)) 

+ E(e^ e > L »l(L G S k , (0, L> < n))E(e^ x(k) ^) 

fee{o,i,2} 

+ ^ E(e f "^h(LeS k ,(6,L)>n)). 

fee{o,i,2} 

Rewriting the left side as 

E(e Me ' L}) l(L e S+)) + J2 E(e M{e ' L) h(L G S k )), 

fee{o,i,2} 

we have 

(1 _ E(e^ 9 ' xi+) »))E(e^ 9 ' L »l(L G 5+, (0, L) < n)) 
fee{o,i,2} 
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Let n go to infinity for this inequality, then the monotone convergence theorem and the 
finiteness of ipi(0) and ^2(6) yield 

o < (i - 7(0) < - VMh) + (lo(0) - iV„(o), 

ke{l,2} 

since G T and / n (x) is nondecreasing in x. The right side of this inequality is finite and 
1 — 7(0) > 0, we must have <p+(0) < 00. 



D Convergence parameter and decay rate 

Let us consider a nonnegative integer valued random variable Z with light tail. Let 

a* = sup{a > 0; Eexp(aZ) < 00}. (D.l) 
a* > by the light tail condition. We give an example such that 

lim x~ l log P(Z > x) = -a* (D.2) 

is not true. One can easily see that 

lim sup — log P(Z > x) = —a*. 

Hence, the problem is the limit infimum. Define the distribution function F of a random 
variable Z by 



F(x) 



1, x < 1, 

e -a*2" 2 n-l < x < 2 n = i 



where F(x) = 1 — F(x), then, for each e > 0, 



lim ^logF(2") = lim ^(-a*2») = -a*, 

lim ^^logF(2" + e) = lim —^—(-a*2 n+1 ) = -2a*. 

n— i>oo 2™ + 6 n— >-oo 2 n + e 



Thus, ( ID.lj) holds, but we have 



liminf — log P(Z > x) = —2a* < —a* = lim sup — log P{Z > x). 

x— >oo X x^-oo X 

The distribution function F(x) is a little tricky because its increasing points are sparse 
as x goes to infinity. However, it is not very difficult to make a small change for it to 
increase at all positive integers. Similar examples are obtained in the literature (e.g., see 
Section 2.3 of |3jJ). 
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E Proof of Lemma 



an 



We first prove 0( A,O °) = r for the three cases separately. For (Dl), suppose that Q^' 00 ^ < 
6 { k k ' c) for k = 1,2. Note that 0j 2,oo) < 0< llOo) and ^ 1,oo) < ef'^ hold by (Dl). Since 

0(l,oo) < 0(A,oo) and e (2,oo) < q(A,oo) ^ at leagt Qne of (l,oo) ^(2,oo) ^ increaged by 

the right hand side of (14. TJ) . This contradicts the supposition. Hence, either ^j 1 ' 00 ^ = #( 1,c ) 
or ^ 2l00) = 6 { 2 2 ' c) holds. Suppose that 0f ,oo) = 0f' c) and ^j, 2 ' ^ < ^' c) . Then, from 
condition (Dl), 6^'°°^ must be increased again by the right hand side of (14. TJ) . This is a 
contradiction. Similarly, it is impossible that #[ 1 ' 00 - ) < 0p c ) but 6^'°°^ = 6^' c \ Thus, we 
must have that 0^'°°^ = Q^'^ for k = 1,2. 

For (D2), we can apply similar arguments as above if we replace B^~ ,c ^ by £,i{6^' c ^)- 
For (D3), we replace of'^ by C, 2 (^'^)- Thus, we get (13.11) for all the three cases since 
r 1 = e ( i L ' oo) andr 2 = e^ oo) . 

The remaining part of this lemma is immediate since <p(0) < oo for all 6 < {0[ A ' n \^ 2 (6[ l 
and for all 6 < (^(O^'™)^^'^) are inductively obtained by Lemma [3.11 



F The proof of Lemma 14.2 



We will use the random walk|l^(f)} introduced in Section [21 We apply the permutation 
arguments in Lemma 5.6 of [5j twice. Then, we have, for any positive integer n and any 



x > 0, 

P(Y(n) ex + A(a), min Y u > 0, min Y 2i > 0\Y(0) = 0) (F.l) 

l<£<n l<£<n 

> 1 -P(Y (n) G x + A (a), min Y 2 £ > 0\Y(0) — 0) 

n \<l<n 

> ^-P(Y(n) ex + A(a)lYYO) = 0). 

n z 

We next note the well known Cramer's theorem (e.g., see Theorem 2 of [4] and Section 
3.5 of 0). 

lim - \ogP{Y{n) G nx + A(o))) = -A(aj), (F.2) 

n->oo n 

where A(x) = sup 0eM2 {(0, x) - log</?(0)}. 

Since the random walk {"K(£)} is identical that of {L(£)} as long as they are inside of 
the quadrant S, (IF.lj) can be written as, for y G S+, 



P{L{n) ex + A(o), do > n\L(0) = y)> —P(Y(n) Ex-y + A(a)\Y(0) = 0). (F.3) 



n 2 



where <Jq = mf{£ > 1; L(£) G dS}. It follows from the representation (14. ip with B = dS 
for the stationary distribution that there are some y G dS and y l G S + such that 
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p(y , Ui) > and, for any m > 1, 
P(L G nc + A(a)) 

= wr^ E E <v)piy> y') E p G nc + A ( a )' ^ > £ \ L ^ = y') 



yedS y'eS^ 



> 



> 



E u (<Tq) 

1 



P(L(m) Enc + A(a), a > m\L{0) = y 1 )v{y )p{y Q , y x ) 
P{Y{m) G nc- Vl + A(o)|y(0) = 0)u(y Q )p(y , yi ) 



m 2 E I/ (a ) 

Thus, for each t > 0, letting m, n — > oo in such a way that n/m — > t, we have 

I 772 1 / 77 

lim — log P(-L G nc + A(a)) > lim log P ( YYm) G m — c — 3/1+ A(a~ 

n-¥oo 77 n->oo r7 777 V 777 

= -^A(tc). 
Since t > can be arbitrary, this implies that 

lim - log P(L G 77C + A(a)) > -inf-A(tc) = - sup{(0, c); 7(0) < 1}, 

71-+00 77 *>0 t 

where the last equality is obtained from Theorem 1 of 0] (see also Theorem 13.5 of [33]). 




Figure 6: The positions of 6, c(cu ), r)(c(cu )) and the cone C max (blue colored area) 
It remains to prove that 6 G" T max implies <p(0) = 00. Define the cone C max as 

C m ax = {XER 2 ;X = S 0^ m&x) + t0( 2 < max ), S ,t>0}. 

If 6 £ C max \ r max , then either 61 > Q^' max ) or 6 2 > 6^' max ^ holds. Hence, <£>(#) = 00 in 
this case by Lemma H~3l Thus, we only need to consider 6 G C max \ r max . 

Since T is a closed convex set that contains 0, there exists a unique 77 G T that 
maximizes (77, c) for each c. Denote this 77 by 77(c). That is, 



(77(c), c) = sup{(T7,c);7(T7) < 1}. 
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Let c(u) = (cosw, sin a;), then the rf(c(u)) continuously moves on <9rnC max from 0( 1 > max ) 
to #( 2 ' max ) as u is increased on (— |7r,7r). Then, as can been seen in Figure El there is an 
ujq G (— \ix, 7r) for G C max \ r max such that Tj(c(uo)) has the same direction as 6, which 
implies that 6 = ar){c{ujQ)) for some a > 1. For this Uq, we have 

(0,c(w o )) = a(r}(c(u )),c(u )) > (r}(c(u )),c(u )). 

Hence, the Markov's inequality, 

(p(0) > e n{e ^ o)) P{L G nc{u Q ) + A(o)), (F.4) 

and (023) yield that 

liminf-log^(0) > (8,c(u )) - (ry(c(w )), c(w )) > 0. 

n— s-oo 7^ 

This concludes ip(6) = oo, which completes the proof. 



G Proof of ( ET2TT ) 



Changing measures for the random walk {Y"(£)} ; we define {Y^ \t)} by 

P(Y (1) (£+1) = (n,j)\Y (1) (£) = (m,i)) = e Tl{n -' m)+T ^- j) P(X {+) = (n-m,i-j)), 
which is well defined since because 7(r) = 1 because of (D2). Note that, for k = 1, 2, 

d 

0=T 



E + 1)| = 0) = E (4 +) e^ (+> )) = ^- 7 (0) 



is finite and positive because of (D2) and the assumption that n < 9i°\ We denote this 
expectation by 

Recall that a[ (n) is the n-th increasing instant of the Markov additive process 
{Z (£)}■ We introduce similar instants for {Y(£)}. Let 

C[ l) (n) = inf{£ > (^(n - 1); Y^{i) - *i (1) (f{ 1} (n - 1)) > 1}. 

For convenience, we also introduce the following events. Let, for 1 < no < n, m > 1 and 
J>0, 



Ki)/~(i) 



7(1)^(1)/ 



"0.' 



inf 

^ 1) (n )<£<^ 1) (n) 

inf r 2 (1) (£)>i,r 2 (1) (c 1 (1) K))=j 
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Since is transient, Z^\n) goes to +00 asn-^oo with probability one. Similarly, for 
k — 1,2, (n) diverges as n — > 00 by fi k > 0. Hence, we can find a positive integer n 
for any i,j ,k ,£ > 1 such that, for any n > n , 

J2 P ^nln(j)\zi 1) (0) = i )>l-e, (G.l) 

j>jo 

E E P (^^') n ft (1) (Cx (1) (n)) = ^}|^ (1) (0) = 0) > 1 - e, (G.2) 

j>jo k>k 

P{Z (l \^\n)) > (£ , 3o )\Z?\o) = i) > 1 - e. (G.3) 
We further choose ri! > n such that, for any n > n 1; 

P(B ^ o (n)|Z«(0)=i)<c (G.4) 
because Z^^^^no)) is finite with probability 1. 
We now compute, for n > n and j > 1, 

= E p ^5» n^ n a)|z (1) (^(n )) = {id)) 



01 



xP(Z (1) (ai 1) K)) = (£,j)|^ 1) (0) = z) 



eT2 A? ,->-» 



= E^ p (^Hn65,„(j)) e _, ; 

xP(Z (1) (ai 1} K)) = V, 3)1^(0) = i), (G.5) 

where the last equality is obtained from the definitions of {Z^(n)} and {Y \n)} and 
the fact that {Y (n)} is a random walk. 

We next consider to apply the renewal theorem. For this let 

^ 1) = ^(n (1) (Ci (1) (i))in (1) (o) = o) 

then is finite because Y^\n) has drifts to +00 by fi^ > and its increments have 
a finite expectation (see Theorem 2.4 in Chapter VIII of [2|). Hence, it follows from the 
renewal theorem that 



lim J2 PpPitP (£)) ' ' 



n 



1=1 MCT 
Obviously, this yields, for each fixed n > 0, 



From (IG.lj) and (]G.4p . it follows that, for sufficiently large n, 

j>jo 
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< e. (G.7) 



Similarly, 



E p fe(n)nc5, n (i)) 



j>jo 



< e. 



(G.8) 



By ( 14.18!) . we can choose sufficiently large jo, k such that 



1 - 



< e, j>jo,k>k 



We then sum both side of ( 1G.5I) for j > j - Further, using ( 1G.2I) and ( 1G.3I) for sufficiently 
large £ , we apply ( 1G.6I) to the sum as n — >■ oo, then letting e -> yield 

lim E P (^no 1 ,n(^) n ^O')!^) =0 = 4)" 



j>jo 



Hence, by (JG\4|l . ( lG~7j) and (IGT8|) . we have 



lim ^P(BiWI?(0) 

j>30 



< e. 



Thus, we have f)4.2ip by letting e | 0. 
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0< L 



maxS 



rnr- 
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